Phase-locking of driven vortex lattices with transverse ac-force and periodic pinning 
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For a vortex lattice moving in a periodic array we show analytically and numerically that a new type 
of phase locking occurs in the presence of a longitudinal dc driving force and a transverse ac driving 
force. This phase locking is distinct from the Shapiro step phase locking found with longitudinal 
ac drives. We show that an increase in critical current and a fundamental phase locked step width 
scale with the square of the driving ac amplitude. Our results should carry over to other systems 
such as vortex motion in Josephson-j unction arrays. 
PACS numbers: 74.60. Ge, 74.60.Jg 



I. INTRODUCTION 

When an external ac drive is applied to a dc driven sys- 
tem interacting with a periodic potential, resonance be- 
tween the external ac frequency and the frequency of mo- 
tion in the periodic potential can give rise to phase lock- 
ing. This phenomenon is found in a wide variety of non- 
linear systems in condensed matter physics. A particu- 
larly well known example is the ac/dc-driven single small 
Josephson-junction and Josephson-junction arrays 
||, where Shapiro steps are observed in the current- 
voltage characteristics. Shapiro type phase locking has 
recently also been studied in experiments Q and theory 
H for vortex motion in two-dimensional superconduc- 
tors with periodic arrays of pinning sites, as well as in 
superconductors where the vortices are driven over a pe- 
riodic potential generated by thickness modulations 
Moreover, driven systems with many degrees of freedom 
in the presence of quenched disorder can exhibit phase- 
locking when there is a dynamically induced periodicity. 
Examples of this are charge density waves |Q and vortex 
lattices in superconductors with random pinning [^|,^). 
In the case of vortices in superconductors with periodic 
pining a vortex can be viewed as an overdamped 

particle moving along a tilted washboard potential where 
the tilt is produced by the dc force from an applied cur- 
rent as well as by the superimposed ac currents in the 
same direction. 

An important difference of vortex arrays from charge 
density wave systems and single degree of freedom sys- 
tems (like a small Josephson junction) is that the dis- 
placement field acting on vortices is a two-dimensional 
vector. This means that displacements can be induced in 
two different directions, and therefore a new kind of phase 
locking is possible when an ac force is applied transverse 
to the direction of the dc force. 

In this work we show analytically and numerically how 
phase locking phenomena may occur when vortices are 
moving in a periodic potential and the ac force is ap- 
plied transverse to the direction of the longitudinally ap- 



plied dc force. The type of phase-locking observed in 
this case is qualitatively different from the Shapiro type 
phase-locking observed when dc and ac forces are in par- 
allel. We first show analytically the possibility of the 
phase locked states for certain commensurate vortex con- 
figurations and predict that the critical current increases 
quadratically with the ac amplitude and that the width of 
some of the phase-locked steps scale as the square of the 
ac amplitude. We also find scaling of the critical current 
as well as the first phase locked region as a function of the 
ac frequency and the pinning geometry. The predictions 
of the perturbation analysis are confirmed with numerical 
simulations. We analyze the validity of the perturbation 
approach and demonstrate how deviations occur. Our 
results suggest that more pronounced transverse phase- 
locking may be observed at lower commensurate mag- 
netic fields, such as B — 1.5£>0 or B = \.2hB§. The 
general model and perturbation results are easily gen- 
eralizable to other systems which exhibit Shapiro step- 
like phase locking such as Josephson-junction arrays and 
Frenkel-Kontorova type models. 

Our studies are directly relevant for several contempo- 
rary efforts since vortex matter interacting with nanos- 
tructured pinning arrays of holes [^0|-jj2| and dots [jl3|jl4| 
has been attracting increasing attention due to the eas- 
ily tunable pinning properties. Pronounced commensu- 
ration effects are observed in these systems when the den- 
sity of vortices matches to integer or fractional multiples 
of the density of pinning sites. In addition to square pin- 
ning arrays, recent experiments have been conducted on 
rectangular pinning arrays |l4f| . Simulations have shown 
several interesting dynamical phases of dc driven vortices 
in periodic pinning systems [ p~5|JT^1 . Imaging [l2| and 
transport experiments along with simulations |l5| ] 

have found that the vortex motion above the first match- 
ing field can occur by the flow of interstitial vortices be- 
tween the pinning sites. These vortices still experience 
a periodic potential created by the vortices located at 
the pinning sites. Recently, phase locking was observed 
for dc and ac driven vortices interacting with periodic 
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pinning at B = 2B r j >1 where is the field at which the 
number of vortices equals the number of pinning sites [Q , 
where the dc and ac drive where in the same direction. 
The values of the voltage response in Ref. j| strongly 
suggest that is is the interstitial vortices that are mobile 
while the vortices on pinning sites remain immobile. 



II. THE MODEL 

We will adopt a simple overdamped particle model for 
describing vortices. The overdamped normalized equa- 
tion of motion for vortex i is 



it = -V i t/({r J -})+f rf , 



(1) 



where r.; = (xi,yi) is the coordinate of vortex i normal- 
ized to the effective 2D penetration depth, A = 2X 2 /d, 
A being the penetration depth and d < A is the sam- 
ple thickness. The energy surface, U , is normalized to 
A v = <5>q/87tA, where 3>o is the flux quantum, and the 
external driving force, = i x z, where i is the nor- 
malized bias current and z is a normalized unit vec- 
tor. Time is normalized to r — aA/A v , where a is the 
Bardeen-Stephen friction coefficient. The energy surface, 
U = U vv ({rj}) + U vp ({rj}), is broken into vortex-vortex 
interactions, U vv ({rj}), and vortex-pinning interactions, 
U vp ({rj}). Vortex-vortex interactions are modeled as 
pair interactions, 



(2) 



i=l j^i 



where m = \vj — r$|, and vortex-pinning interactions are 
modeled by, 



i=l j=l 



(3) 



where N v and N p are the number of simulated vortices 
and pinning centers, respectively. 

We are performing simulations with periodic boundary 
conditions (periodicity L x , L y ), and with a logarithmic 
interaction potential (u vv (r^ ) = — In ry ) valid when A 
is of the order of the sample size. 

Combining the logarithmic interactions with the peri- 
odic boundary conditions, result in Jl7| 



,,,.„(,,,.,,,,) ' ( 7^1 ~ 7T ~r~ (f- 



(4) 
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g cosh ( 2 ^(^+n)) 
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cosh 



where (acy, are the cartesian coordinates, r?- = xf, + 
yfj, and C(-) is a function determined by the energy nor- 
malization of the long range interactions. 



We have previously demonstrated that phase-locking 
of the vortex motion can exist if a certain periodicity of 
the global pinning potential is present. This was shown 
for driving forces of the form, = (77 + 6 sin f2i)x, where 
x is a unit vector in the x-direction. Analogous to the 
well-known Shapiro steps Q| in the current-voltage char- 
acteristics of an ac/dc-driven Josephson junction, it was 
demonstrated that nonzero intervals of rj will result in 
phase-locking of vortex motion with vortex velocities res- 
onating with the external ac-field. We will in this paper 
demonstrate that one can also expect phase-locking in 
the transversely ac-driven case, i.e., for 



e sin ilt , 



(5) 



where rj and e are the magnitudes of dc and ac forces and 
f2 is the frequency of the ac drive. The relevant equations 
of motion arc therefore, 



±i + i i u{{T]]) = ri 

d 



III. ANALYTICAL RESULTS 



(6) 
(7) 



We will consider a system with periodic pinning in a 
rectangular periodic lattice with plaquette area A = l x l v , 
where l x and l y are the dimensions of the rectangular pla- 
quettes. We will assume that the applied magnetic field 
is larger than the first matching field, such that all pin- 
ning centers are occupied by trapped vortices. We will 
in the following assume that the trapped vortices remain 
trapped at the pinning sites and immobile at all times. 
Thus the interstitial vortices move in a strictly periodic 
lattice of immobile vortices with periodicity (l x , l y ). We 
can then perform a perturbational analysis of the inter- 
stitial vortex behavior through the following assumption: 
any interstitial vortices moves such that the forces from 
other interstitial vortices cancel by symmetry B. Thus, 
analyzing the behavior of the interstitial (mobile) vor- 
tices, we will adopt the potential function for a single 
interstitial vortex moving in a lattice of pinned vortices 
(at positions (nl x , (m + as given in Eq. (|^). Thus, 

the interstitial vortex position is subject to a potential 
as in Eq. (H) where the coordinate is the vortex posi- 
tion, (Xi,yi), and with the lattice constants here being 
the distances between pinning centers, (l x , l y ), instead of 
the computational system dimensions, (L x ,L y ). 



A. Critical dc-force without ac-drive 



The first result to evaluate for the purely dc-driven 
case is the "critical dc current", r]i°\ below which the 
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interstitial vortex is trapped by the pinned vortex lat- 
tice. As 77 is increased from 7/ = 0, the interstitial vortex 
will find an equilibrium point for which the pinned vortex 
lattice provides a cancelling force to the dc bias. How- 
ever, for a given (critical) value of 77, namely rjc°\ the 
pinned vortex lattice can no longer provide enough force 
to resist motion of the interstitial, which therefore be- 
gins to propagate in the x-direction. Assuming that the 
pinned vortices do not move and that an interstitial vor- 
tex only interacts with pinned vortices, this critical force 
can be found as the maximum gradient value of Eq. (^) 
for yij = mly. The corresponding equation of motion is 
given by Eq. (^J): 



Xi + — ^ 



sin 



cosh ( 2 ^(i+n)) - cos (2^) 



(8) 



This equation can produce a critical current by requir- 
ing in = and optimizing the left hand side through the 
varying the position, a;,-. In general, this must be done nu- 
merically. The result is shown in Fig. (1), where the crit- 
ical dc force (thick solid) is shown together with the op- 
timized value of the equilibrium value of Xi = . In the 
limit of sech{irl y /l x ) <C 1, we can derive the approximate 
expression of the critical current by noticing that only 
two terms (n — —1,0) are required to describe the criti- 
cal current, which is then always given by X{ — x* = jl x : 



= /cSCCh ( 7T j- 



, 7T l x 



(9) 
(10) 



The approximate expression, Eq. (|9j), is shown in figure 
1 as a thin solid curve, validating that the interaction 
potential between interstitial and pinned vortices can be 
simplified to the terms n = —1,0 for l x ^1.5l y . One can 
also see that in the extreme opposite limit, l x 3> l y , the 
critical dc force, r]i°\ can be approximated by tt/I v . This 
is shown as a dashed line in figure 1. 

It is here important to remember that the analysis, 
which is based on a magnetic field of "matching field 
plus one flux quantum" , is actually valid for magnetic 
fields, < B < 2f?0, for which the vortex configuration 
results on forces between interstitial vortices cancelling 
due to symmetry; such magnetic fields can be, e.g., B = 
2I?0, jB^, ■ ■ ■. 



B. Transverse ac-drive included 



Let us assume that sech 



(j r T L ^J C 1, so that we can 
write the approximate equations of motion for the inter- 
stitial vortices as, 



(0) 



;k 2 yi 



sin (xik) = r/ 



27r • r, 

In" I? , 



(11) 

(12) 



where we have retained terms in yi only to lowest order. 
We have required 7r|j/i <C l x in order to obtain Eqs. ( pd| ) 
and @. 

The solution to the second of those equations is easily 
found, 



e sin(fit + 9) 



with 



Wl 



<!/,/, 



tan 



2tt 

£ll X ly 



Inserting Eq. (O) into Eq. (fil) yields, 



(0) 



1 + 

4 



sin kx; 



(13) 
(14) 
(15) 

(16) 



7/<; 0) -fc 2 £ 2 cos2(fi7j + 9) sin fca 



>1 



This equation is the effective equation of motion for the 
vortex behavior in the longitudinal (x) direction given 
an ac-force in the transverse (y) direction, derived within 
the approximations listed above. The equation is equiv- 
alent to that of an overdamped pendulum (kxi being the 
pendulum phase) with dc-torque (77) and a pivot verti- 
cally oscillating with frequency 2il and amplitude pro- 
portional to e 2 - i.e., an overdamped equivalent of the 
classic Kapitza problem p"8| ] for underdamped and para- 
metrically driven pendula. 

We will now consider a few separate cases of vortex 
responses to the transverse ac-force: 
(x.i) = <^> Xi(t) — x : We will omit the last term on the 
left hand side of Eq. (|l6|). The remaining (dc) terms are: 

77c sin kxo =77 (17) 

8r, c (18) 



2tt 



1 



1 + 7 fc2 ^ 
4 
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\n) ft 



1 
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\QUyJ 



W(S1 



2 

x 1 



2tt 

ill, I, 



(19) 



where rj c is the critical dc force. It is here worthwhile to 
notice that the critical current, ?7 C , increases (Sr/ C posi- 
tive) quadratically with the transverse ac-amplitude, e. 
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This is in direct contrast to the case when the ac-drive is 
longitudinal, in which case the critical dc force decreases 
quadratically 0. We remember that validity of the ex- 
pressions requires 7re <C Cll x . 

(ii) > <^ 1 77 1 > r) c : Without the parametric term in 
Eq. (|l6|), the solution, for rj 2 > r/ 2 , is given by ||] 



''0 



Vc > 



1 - 



1 — 2a C os kvnt ' 

1)0 u 



(20) 
(21) 



where rjo is the dc-force when the parametric term in 
Eq. (16) is neglected. Let us assume the following simple 
ansatz for the vortex motion, 



Xi(t) = v Q t + Xq 



(22) 



where xq is a constant. Inserting this ansatz into the 
effective equation of motion, Eq. (|l6|), while only main- 
taining dc terms, yields, 

v o — 5r] c (cos2(ilt + 9) sin kxi(t)) = rj = 770 + Srj . (23) 



With the ansatz of Eq. (|22|) we can therefore expect the 
term (•) to contribute to 8f] if kvo = 2f2 (see figure 2), 
and the resulting relationship between the internal phase, 
kxo — 29, and the dc current, rjo + Srj, is: 



20,/k — —6rj c sm(kxo 



29) = 770 + St) 



(24) 



As the phase, kxo — 29, can be adjusted to balance this 
equation for different choices of £77, we can argue for 
a nonzero range, Arj, in dc-force for which the aver- 
age speed of the interstitial vortices is unchanged. This 
phase-locking range has the magnitude: 



Arj = Sf] c . 



(25) 



Thus, we can predict phase-locking in the transversely 
ac-driven vortex system and the predicted total range in 
phase-locking is equal to the increase in the critical de- 
force due to the ac-drive. This prediction is correct up 
to and including terms cx e 2 . 

This range in bias current for which the average speed 
(voltage) of the interstitial vortex is constant will man- 
ifest itself as a step in the dc current-voltage character- 
istics of the system. It is important to emphasize that 
this step is very different in origin from the Shapiro steps 
H recently demonstrated for longitudinal ac-drive. The 
essential difference lies in the fact that the Shapiro steps 
arise from a "direct" driving term in an effective pendu- 
lum equation, whereas the present phase-locking for the 
transverse ac-drive arises from an effective parametric ac- 
driving term in the longitudinal equation of motion. 

If one chooses a better (and more complicated) ansatz 
for the vortex motion (e.g., Eq. ( ^T| ) instead of Eq. (|2^)) it 



becomes evident that an e 2 phase-locked step also exists 
for kvo — £1 and that higher order phase-locked steps 
may exist at any sub and super harmonic of the driving 
frequency. However, we will not go into detail with other 
phase-locked modes in this presentation. 



IV. NUMERICAL SIMULATIONS 

In order to validate the analysis of Sec. Ill predicting 
some basic effects of a transverse ac- force on vortices, we 
have conducted numerical simulations of driven intersti- 
tial vortices in a rectangular lattice of pinned vortices. 

Figure 3 shows a series of simulated (normalized) cur- 
rent (dc-force) voltage (average vortex velocity, vq) for 
different values of the transverse ac-drive. The system is 
a square array of 4 x 4 pinned vortices with periodicity 
(l x , l y )=(2, 2) in a computational simulation box of size 
(L x ,L y ) = 4l x ,Al y ). The simulated magnetic field corre- 
sponds to B/B<f, = 1.5 such that each pinning center is 
occupied by a vortex and every other plaquette has an 
interstitial vortex. The applied normalized frequency is 
= 4. The four curves (shifted vertically for clarity) in 
figure 3 are for (from top) e = 0, 1, 2, 3. According to the 
analysis above, the critical dc-force, rj c , resulting in onset 
of voltage (vortex transport) should increase quadrati- 
cally for increasing e. This is clearly visible from figure 
3. A phase-locked step in dc bias current should develop 
for increasing e around kvo = 2Q. This is also clearly 
visible for rj m 2.5. As the ac amplitude is increased we 
observe more steps in the IV characteristics. However, 
while the simplest of these, kvo = f2, can be analyzed in 
some detail, we will be concerned with only the steps at 
kvg = and kvo — 2il here. 

Two characteristically different vortex evolutions are 
shown as insets to figure 3. At bottom right we show a 
snapshot (•) of a vortex configuration within the param- 
eter range of the phase-locked step for e = 2 together 
with a long time trace (thin line) of the vortex trajecto- 
ries. The vortex paths are obviously periodic and the in- 
terstitial vortices seem to move in a geometrically simple 
configuration where all interactions between the intersti- 
tials cancel due to symmetry. The top left inset shows 
the traces outside of the phase-locked region. Here we 
observe the pinned vortices as static dots, while the inter- 
stitials move in a chaotic or quasiperiodic band between 
the pinning sites. 

A detailed investigation of the validity of the above 
perturbation results was performed for several different 
sets of system parameters. For a system of a single inter- 
stitial vortex, we show in figure 4 the detailed compar- 
isons between the analytical predictions (lines), Eq. jis|), 
and numerical simulations (markers) of the increase, dr] c , 
in critical dc force as a function of ac-amplitude. Fig- 
ure 4a shows data for square pinning arrays of different 
plaquette areas (l x ly = 4, 16) and with different driving 
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frequencies (O = 4,8), while figure 4b shows compar- 
isons for rectangular systems (l x /l y = 5,2) of one pla- 
quette area (l x l y = 4) for different frequencies (O = 4, 8). 
The details of the parameter sets are given in the cap- 
tion. The agreement between simulation and analysis 
is obviously overall good for these different parameters. 
Certainly the quadratic relationship between the increase 
of the critical current and the ac-amplitude is confirmed. 
The quantitative agreement is also rather good. The sim- 
ulation results seem underestimated by the predicted val- 
ues. However, this is to be expected since the analysis 
only includes the transverse vortex motion as moving in 
a harmonic potential. The vortex model will provide a 
stronger nonlinearity for larger ac-amplitudes, resulting 
in a larger than predicted increase of the critical current. 

For the same parameter sets as shown in figure 4, we 
have performed simulations of the phase-locked step at 
kvo = 2Q and compared the results to the predicted 
range, A77, given in Eq. (p5|). The results, shown in figure 
5, are also here in good agreement with the predictions. 
Thus, we are confident that the above analysis of the 
single interstitial behavior can be a useful tool for under- 
standing phase-locking in this system. 

Figure 6 shows, again for the same parameter sets as in 
figure 4, the comparison between the magnitude of the 
phase-locked step, A?;, and the predicted value, but in 
this case for B — l.hB^ i.e., for one interstitial vortex 
per two plaquettes of pinned vortices. Once again we 
find rather good agreement between simulations and the 
simple prediction. However, we notice that rectangular 
systems with l x — 2l y (top curves in figure 6b) cannot 
produce a measurable phase-locked step for very small 
ac-amplitudes. 

If one simulates the system for the same parameter 
values, but for B — 2B^ (not shown), the agreement will 
get rather poor for all rectangular cases where l x = 2l y , 
and even for some of the square array cases as well. The 
reason for these discrepancies is of course found in the 
assumptions behind the perturbation analysis. Partic- 
ularly, the validity of the assumption of noninteracting 
interstitial vortices will break down for densely popu- 
lated (B = 2B^) systems since the vortices are mutu- 
ally repulsive. The result is sketched in figure 7, where 
the vortex repulsion is deforming the geometrically sim- 
ple lattice into a configuration where the effective inter- 
stitial vortex-vortex interaction is non-zero. The lat- 
tice deformation is characterized by the two parame- 
ters, £( x > and Fi gure 8 shows simulation results 



for L x = Al x = L y 



= = 



8, Q = 4, and B = 2B d 



The circular markers (•) show the simulated range in 
phase-locking and the solid straight line represents the 
prediction of Eq. (p5|). It is apparent that the agreement 
between simulations and prediction is only fair (A77 oc e 2 ) 
for relatively large values of e. Measuring the average (in 
time and space ) "phase difference", (\£ {x) \) and 
we can correlate the magnitude of lattice deformation 



with the comparison between simulation and analysis of 
phase-locking. The measures of lattice deformation are 
defined as, 



O) 

i 

Av) 



(xi mod l x ) 

(UiVHOAly) 



,(x) 
>i,3 



Jx) _ Jx) 
£7 M 



It. 



Jx) _ Jx) 





>-< 








Ay) _ Av) 



Ay) _ Ay) 



A x ) _ A x ) 
Av) Ay) 



> 1 



Ay) _ Av) 

^7 



< i3L 

— 2 



> 



^ 2 





A?) 


2 )«-( 






Av)\ 







(26) 
(27) 

(28) 



(29) 

(30) 
(31) 

(32) 

(33) 
(34) 



Finally, these averages are averaged over the phase-locked 
step to provide a single measure for a given e ((\^ x ' v ^\) 
and a x , y are typically smallest near the center of the lock- 
ing range). Thus, < (l^^l) < l x , y /2 is a measure of 
the average spatial lattice deformation and a XiV is a mea- 
sure of the average temporal fluctuations in the lattice 



deformation, 



i,3 



The minimum averaged lattice deformation parame- 
ters for a given phase-locked step are shown in figure 8 as 
solid markers (□: A: {\(,^ |)). We observe large 

lattice deformations at small ac-amplitudes, where the 
agreement between the observed range of phase-locking 
is poor, while the deformations "collapse" for larger val- 
ues of e, thereby giving rise to better agreement between 
simulations and analysis. We have also shown the average 
of the temporal fluctuations, a x and a y , in order to illus- 
trate that the lattice deformations depend on time when 
<|£^|) and (\^ x) \) are non-zero. Thus, it indicates that 
internal modes of the lattice deformation are excited and 
may contribute significantly to the phase-locking range in 
the low ac-amplitude cases. For large ac-amplitudes we 
observe a saturation of the range in phase-locking accom- 
panied by reappearance of the deformation lattice. This 
is consistent with chaotic behavior which are expected 
for strong ac-driven nonlinear systems. 
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V. CONCLUSION 

We have developed a simple analysis of the behav- 
ior of interstitial vortices in systems with periodic pin- 
ning and transverse ac-drive. The result indicates that 
the dc-driven longitudinal motion can phase-lock to the 
transverse ac-signal and that the range of phase-locking 
(in dc-current) is quadratic in the ac-amplitude and we 
have developed a quantitative expression providing de- 
tailed dependencies of also other relevant system param- 
eters, such as pinning geometry and driving frequency. 
We have further demonstrated that the critical current 
increases quadratically with the ac-amplitude. The per- 
turbation results have been validated by numerical sim- 
ulations which show good agreement with the analyt- 
ical predictions in the expected range of system pa- 
rameters. The mechanism of phase-locking discussed in 
this paper is distinct from phase-locking to a longitu- 
dinally applied ac-signal. The latter case was studied 
in U and exhibits qualitatively different responses to 
ac-perturbations, such as decreasing critical dc-current 
with increasing ac-amplitude and phase-locked steps that 
grow linearly with ac-amplitudc. The results shown in 
this paper are for commensurate fields. Incommensurate 
fields, where no simple geometrical relationship can exist 
between the pinned and interstitial lattices, are charac- 
terized by non-cancelling interactions between intersti- 
tial vortices. Thus, phase-locking ranges for the non- 
commensurate fields usually have magnitudes less than 
predicted by the single interstitial analysis. 

Our results and analyses indicate that the most impor- 
tant phase-locking appears at relatively small magnetic 
fields (with respects to the first matching field) where 
the inter-vortex repulsion is not deforming the intersti- 
tial vortex lattice. 

The predictions in this paper should be directly ap- 
plicable for experimental verification in superconductors 
(or Josephson junction arrays) where dc and ac fields are 
orthogonal. 
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FIG. 6. Magnitude, Arj, of the phase-locked step at 
kvo — 20 as a function of ac-amplitude, e. Markers repre- 
sent numerical simulations of equations (|^) and ([? ) . Solid 
lines are the corresponding predictions from Eq. (125). Simu- 
lations are conducted for B = 1.5£>0, L 
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FIG. 1. Static critical force, ??c x, for e = as a function of 
rectangular periodic pinning aspect ratio. Thick lines repre- 
sent the maximum force derived from Eq. (^) for the optimal 
position, x„, along the x-axis. Thin line represent Eq. (|^), 
which is valid for sech(7rZ H /l x ) <C 1, dashed line represent 
n/ly, which is asymptotically correct for l x 3> l y - 

FIG. 2. Sketch of interstitial vortex (•) trajectory in rect- 
angular lattice of pinned vortices (o). Upper part of sketch 
shows a case where (xi) = while the lower part shows a 
case where (±i) = l x Q/ir. Ac-force is transverse, dc-force is 
longitudinal. 

FIG. 3. Simulated IV characteristics for different values 
of transverse ac-amplitude, e = 0,1,2,3 (top down). IV 
curves are vertically offset for clarity. System parameters are: 
L x = L y = 4l x = 4l y — 8, O = 4, and B = 1.5-B^,. Lower 
inset shows vortex trajectories in a phase-locked state while 
upper inset shows trajectories outside a phase locked regime. 



FIG. 7. Sketch of interstitial vortex lattice deformed due 
to the repulsive vortex- vortex interaction at B = 2B^,. The 
parameters, and ^ y ', are measures of the magnitude of 
deformation from a rectangular lattice. 

FIG. 8. Magnitude, Arj, of the phase-locked step 
at kvo — 2Q. as a function of ac-amplitude, e, for 
L x = 4l x = L y = il y = 8, fl = 4, and b = 2Bj, (•). 
Markers represent numerical simulations of equations M) and 
(0). Solid line represent the corresponding prediction from 
Eq. (Ba). Also shown are the lattice deformation measures, 
(l^^'J) and (|C' y '|), as well as their standard deviations, a x 
and a y . 



FIG. 4. Increase of critical dc current, 8r] c as a function 
of ac-amplitude, e. Markers represent numerical simulations 



of equations (|6|) and (7J). 
predictions from Eq. (191 
single interstitial vortex. 
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4 and ft = 8. 



(a) A: l x 
8; □: L = 



I y 



l X ly 

O: l x = 2l y = 4 and Q = 8 
2l x =l y =4 and Q = 8. 



(b) A: l x = 
□ : 2l x = I 



2ly 



= 2 and = 4 
4 and Q = 4; o 
4 and = 4 
4 and Q, = 4; o 



FIG. 5. Magnitude, A77, of the phase-locked step at 
kvo — 20. as a function of ac-amplitude, e. Markers represent 
numerical simulations of equations ([]) and (Q) . Solid lines are 
the corresponding predictions from Eq. (^H|). Simulations are 
conducted for a single interstitial vortex, (a) A: l x — l y = 2 
and Q. = 4; O: ^ = ly = 2 and fi = 8; □: Zz = i s = 4 and 
Q, = 4; o: Z* = Z y = 4 and = 8. (b) A: Z x = 2l y = 4 and 
= 4; O: Z,; = 2Z y = 4 and Q = 8; □: 2Z,; = l y = 4 and 
= 4; o: 2ZZ = l y = 4 and = 8. 
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